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A CALCULATION  STUDY  OF  WING-AILERON  FLUTTER  IN 
TWO  DEGREES  OF  FREEDOM  FOR  TWO-DIMENSIONAL 
SUPERSONIC  FLOW 

By  Donald  S.  Woolston  and  Vera  Huckel 
SUMMARY 


Results  are  presented  of  some  sample  calculations  made  for  the 
bending-aileron  and  torsion-aileron  flutter  of  an  idealized  wing-aileron 
system  (infinite  aspect  ratio,  full-span  aileron  hinged  at  its  leading 
edge)  for  several  values  of  Mach  number  in  supersonic  flow  and  also  for 
some  subsonic  Mach  numbers.  It  is  emphasized  that  the  results  are  sub- 
ject to  various  limitations  and  are  not  intended  to  be  applicable  to 
particular  configurations.  They  may  serve,  however,  to  provide  pre- 
liminary knowledge  of  the  influence  of  variations  in  certain  parameters 
and  to  indicate  some  of  the  differences  in  the  trends  exhibited  by  the 
calculations  for  high  speeds  and  those  trends  known  to  exist  for  low 
speeds. 

The  more  influential  parameters  appear  to  be:  mass  balance  of  the 

control  surface,  control-surface  frequency  relative  to  wing  frequency, 
and  structural  damping.  An  interesting  result  is  the  indication  of  an 
apparent  reverse  effect  of  mass  balance  at  the  transonic  speeds  for  some 
wing-aileron  configurations.  A high  value  of  control-surface  frequency 
appears  beneficial  as  does  structurrl  damping  of  the  system.  In  several 
of  the  calculations,  the  possibility  of  flutter  in  a single  degree  of 
freedom  is  indicated. 


INTRODUCTION 


Numerous  investigations,  both  analytical  and  experimental,  of  the 
problem  of  flutter  involving  control  surfaces  have  been  made  in  order 
to  determine  what  parameters  influence  this  type  of  flutter  and  to  find 
means  of  preventing  it.  These  investigations  have  been  concerned  chiefly 
with  flutter  at  low  speeds.  The  present  paper  aims  at  contributing  to 
the  subject  through  calculations  made  by  analytical  consideration  of  the 
supersonic  speed  range.  Since  the  analysis  has  been  idealized  and  has 
many  limitations,  sane  remarks  emphasizing  the  nature  of  the  limitations 
follow  this  introduction  and  the  list  of  symbols. 
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Some  of  the  factors  affecting  control-surface  flutter  at  low  speeds 
have  long  been  known.  Even  before  the  development  of  the  field  of 
unsteady  aerodynamics,  some  conclusions  had  been  reached  by  Von  Baumhauer 
and  Koning  (ref.  l)  on  the  influence  of  control-surface  mass  balance. 
Similar  conclusions  were  reached  by  Frazer  and  Duncan  in  reference  2. 
later  investigations  of  control -surface  flutter  include  the  treatment  of 
Theodorsen  and  Garrick  (ref.  3) , where  the  results  are  presented  of  a 
systematic  analytical  survey  of  the  effects  of  various  parameters  on 
the  flutter  characteristics  of  a wing-aileron  system  in  incompressible 
flow.  Other  investigations  have  been  made  by  Van  de  Vooren  (ref.  4), 
numerous  British  investigators  (see,  for  example,  ref.  5),  Dorr  (ref,  6), 
and  others.  This  work  has  indicated  the  critical  dependence  of  control- 
surface  flutter  on  such  parameters  as  control- surface  unbalance,  struc- 
tural damping,  and  the  relative  stiffnesses  in  wing  bending,  wing  torsion, 
and  aileron  motion. 

The  criteria  determined  analytically  for  the  case  of  incompressible 
flow  have  generally  been  carried  over,  together  with  experience,  to  form 
the  basis  for  criteria  for  the  prevention  of  control-surface  flutter  in 
the  subsonic  case.  Although  subsonic  aerodynamic  coefficients  relating 
to  the  control  surface  have  been  available  in,  for  example,  reference  7, 
in  certain  British  sources,  and  more  recently  in  more  extensive  form  in 
reference  8,  no  systematic  calculations  using  these  coefficients  have 
been  published. 

With  the  attainment  of  supersonic  speeds  in  flight  the  question  of 
the  applicability  of  standards  set  for  the  subsonic  case  to  the  super- 
sonic case  has  become  increasingly  important.  Little  practical  experi- 
mental knowledge  exists  and  until  recently  no  theoretical  aerodynamic 
coefficients  relating  to  control-surface  flutter  at  supersonic  speeds 
were  available.  These  coefficients  have  now  been  tabulated  for  a wing 
in  two-dimensional  supersonic  flow  in  reference  9 as  an  extension  of 
reference  10.  Coefficients  for  the  sonic  case  have  been  developed  and 
tabulated  in  reference  11. 

In  the  present  paper  the  aerodynamic  coefficients  of  references  9, 

10,  and  11  have  been  employed  with  the  analytical  approach  of  reference  10 
in  the  performance  of  some  sample  calculations  for  wing-aileron  flutter 
in  two  degrees  of  freedom  at  supersonic  speeds.  For  comparison,  a few 
examples  for  subsonic  flow  are  included.  The  investigation  applies  to 
an  infinite-aspect-ratio  wing-aileron  system  with  an  aileron  hinged  at 
its  leading  edge  and  is  not  concerned  with  any  existing  configuration. 
Control-surface  configurations  for  supersonic  flight  are  not  yet  stand- 
ardized and  hence  the  investigation  has  been  limited  to  an  idealized 
system  with  the  aim  of  gaining  some  indications  of  what  to  look  out  for 
and  what  to  expect  in  control-surface  flutter  at  supersonic  speeds  as 
compared  with  what  is  known  for  the  subsonic  case.  It  is  realized  that 
similar  calculations  may  exist  elsewhere  in  unpublished  form  and,  further, 
that  a solution  to  the  complete  problem  will  require  the  consideration 
of  many  factors  not  included  in  the  present  work. 
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SYMBOLS 


b 

c 

h 


I 


P 


M 

Mn»  % 
(n  = 1,2,.. 


rP 


one -half  chord  of  wing-aileron  system 

velocity  of  sound 

structural  damping  coefficients 

vertical  displacement  of  axis  of  rotation 

moment  of  inertia  of  wing-aileron  combination  about  elastic 
axis  per  unit  span  length 

moment  of  inertia  of  aileron  about  x^  (hinge  axis)  per 
unit  span  length 

mass  of  wing-aileron  combination  per  unit  span  length 
Mach  number 

aerodynamic  coefficients  of  lift  and  moment  (see 
refs.  9 and  10) 

radius  of  gyration  of  wing-aileron  combination  referred 
to  x0,  yio/nib2 

reduced  radius  of  gyration  of  aileron  referred  to  x^, 

^ Ip/mb2 


Sa  static  moment  of  wing-aileron  combination  per  unit  span 

length  referred  to  xQ 

Sn  static  monent  of  aileron  per  unit  span  length  referred 

to  x-^ 

v speed  of  forward  motion  (or  velocity  of  main  stream) 

xQ  coordinate  of  elastic  axis  measured  from  leading  edge  of 

wing,  referred  to  chord  2b  as  reference  length 

coordinate  of  aileron  hinge  measured  from  leading  edge  of 
wing,  referred  to  chord  2b  as  reference  length 

*a  location  of  center  of  gravity  of  wing-aileron  combination 

referred  to  xQ 


f 


l 

i 


-4 


r! 


xu 

a 

0 

H 

0) 

<% 


0), 


3 


NACA  TN  3160 

reduced  location  of  center  of  gravity  of  aileron  referred 
to  x^,  Sp/mb 

angular  displacement  about  rotation  point 
angular  displacement  of  aileron  about  aileron  hinge 
mass  ratio  parameter,  m/kpb 2 
circular  frequency  at  flutter 

natural  circular  frequency  of  wing  in  vertical  translation 

natural  circular  frequency  of  torsional  vibrations  about 
elastic  axis 

natural  circular  frequency  of  torsional  vibrations  of 
aileron  about 

mass  density  of  air 


fihx  = ^(a)hH2(1  + ^h) 

fiax  = ^ra2(%/“)2(1  + ^ 
fipx  = nr32^a)3/a))2(l  + igp) 


REMARKS  EMPHASIZING  LIMITATIONS 


The  calculations  discussed  in  this  report  were  undertaken  in  order 
to  gain  sane  insight  into  the  problem  of  flutter  involving  control  sur- 
faces at  supersonic  speeds,  a field  in  which  very  few  published  results 
exist,  and  were  intended  in  part  to  serve  as  a possible  guide  to  experi- 
mental investigations  of  the  problem.  It  is  felt  that  the  results  of 
the  numerical  calculations  mJ ght  be  of  interest  to  other  investigators. 

It  should  be  clearly  kept  in  mind  that  the  results  presented  are  subject 
to  various  limitations  and  are  not  intended  for  application  to  an  actual 
configuration.  As  indicated  in  the  introduction,  the  calculations  were 
performed  for  the  case  of  two-dimensional  flow  because  of  the  availability 
of  the  coefficients  and  were  applied  to  a highly  idealized  wing -aileron 
system  (infinite  aspect  ratio,  full-span  aileron  hinged  at  its  leading 
edge).  Parameter  variations  treated  are  rather  limited.  The  results 
may  be  considerably  altered  by  such  factors  as  finite  span,  thickness, 
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additional  degrees  of  freedom,  and  other  parameters  which  were  not  con- 
sidered. The  results  should  serve,  however,  to  provide  a preliminary 
knowledge  of  the  influence  of  variations  in  certain  parameters  and  to 
Indicate  some  of  the  differences  in  trends  exhibited  by  the  calculations 
for  high  speeds  and  the  trends  known  to  exist  for  low  speeds. 


SYSTEM  CONSIDERED  AND  METHOD  OP  ANALYSIS 


Consider  a wing-aileron  system  having  the  degrees  of  freedom  of 
vertical  translation  h,  angular  displacement  a,  and  aileron  rotation  0. 
This  system  will  be  replaced  in  the  analytical  approach  employed  herein 
by  the  following  system  of  infinite  aspect  ratio 


where  the  chard  and  other  dimensions  are  taken  as  those  of  a repre- 
sentative wing  station.  The  stiffnesses  of  the  linear  spring  and  the 
two  torsion  springs  are  chosen  to  give  natural  frequencies  which  are 
the  same  as  the  bending,  torsion,  and  aileron  natural  frequencies  of 
the  wing  system. 

It  is  shown  in  reference  10  that  for  a system  of  this  type  the 
flutter  determinant  is 
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where 


Ach  = 

M - p + 

+ iLg 

^ca  = 

+ ^ 

+ ity 

Acp  = 

+ I5 

+ 1^6 

Aah  = 

-nxa  + M1 

+ iMg 

^a a 

< - W*  + Mj 

+ iM4 

= 

r32  - 2(xl  " 

Xo)x^] 

+ ^ + 

iMg 

<? 

& 

(1 

-HXp  + Nx 

+ iN2 

Aba 

-4 

Es2^ 

>(xi  - 

xo)x0_ 

+ Nj  + 

i% 

Abp  = 

<V 

< - V*ff 

- + n5 

+ iN6 

The  are  functions  of  the  structural  properties  of  the  wing 

and  the  flutter  frequency  cu  through  the  relations 

nhx  = ^K/co)2(1  + ^h) 

naX  = ^ra2(ak/a>)2(1  + laa) 

Sx  = ^r02(a>p/“)2(1  + ige) 
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lhe  L's,  M's,  and  N's  are  aerodynamic  force  and  moment  coefficients 
defined  and  tabulated  in  reference  9 and  10  for  the  supersonic  case,  and 
are  given  for  the  sonic  case  in  reference  11.  The  coefficients  for  sub- 
sonic compressible  flow  can  be  obtained  for  certain  limited  cases  in 
references  7 and  8 and  for  incompressible  flow  from  reference  12.  A 
discussion  of  the  significance  of  the  structural  parameters  and  of  the 
numerical  values  assigned  to  them  is  given  in  a subsequent  section. 

The  calculations  to  follow  have  been  limited  to  the  two  subcases 
which  involve  two  degrees  of  freedom,  that  is,  flutter  involving  wing 
bending  and  aileron  rotation,  and  flutter  involving  wing  torsion  and 
aileron  rotation.  It  is  recognized  that,  in  cases  of  high  coupling 
between  bending  and  torsion,  a three -degree -of -freedom  analysis  would 
probably  be  required  to  give  a complete  picture  of  the  situation.  In 
cases  which  are  not  highly  coupled,  however,  consideration  of  the  sub- 
cases should  lead  to  satisfactory  results,  at  least  with  regard  to 
trends . 


DISCUSSION  OF  PARAMETERS  FOR  THE  CASES  TREATED 


The  various  parameters  to  be  considered  will  be  Mach  number  M, 
mass  ratio  parameter  (i,  aileron-hinge  location  x^,  aileron  radius  of 

gyration  r^,  unbalance  parameter  x^,  frequency  ratios  Jcu^  and 

GDp/o^,  wing  radius  of  gyration  ra2,  elastic-axis  location  xQ,  and 

structural  damping  gh,  g^.  A table  of  values  of  these  parameters 
employed  in  the  calculations  is  presented  as  table  I. 

A variation  in  the  mass  ratio  parameter  p.  may  be  considered  as  a 
variaton  either  in  altitude  or  in  wing  mass.  The  values  of  u = 15.708 
and  157*08  have  been  employed  in  the  calculations  and,  for  a particular 
configuration,  might  represent  altitudes  of  sea  level  and  approximately 
50,000  feet,  respectively. 

The  parameter  x^  is  the  coordinate  of  the  aileron  hinge  relative 
to  the  leading  edge  of  the  wing.  A value  of  x^  = 0.8,  for  the  hinge 

at  80  percent  chord,  has  been  used  in  the  majority  of  the  calculations. 
Some  results  are  given  for  the  hinge  at  50  percent  chord,  or  x^  = 0.5. 

The  parameter  r^2  as  defined  is  related  to  the  radius  of  gyration 

of  the  aileron  referred  to  the  hinge  and  depends  in  part  on  the  ratio  of 
aileron  mass  to  wing  mass  and  on  the  aileron  mass  distribution.  For  the 
case  where  x^  = 0.8,  if  a ratio  of  aileron  mass  per  unit  length  to  wing 
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mass  per  unit  length  of  1/10  is  assumed,  for  a uniform  mass  distribu- 
tion, r p2  would  he  approximately  0.005*  whereas,  for  a triangular  mass 

distribution,  rp2  would  be  approximately  0.003.  An  average  value  of 

rp2  = 0.004  has  been  used  in  most  of  the  calculations.  For  the  case 
where  x^  = 0.5  the  assumption  of  a mass  ratio  of  1/3  would  lead  to  a 
similar  average  value  for  rp2  of  about  1/10. 

The  parameter  Xp  is  a reduced  parameter  which  also  involves  the 

mass  ratio  and  is  related  to  the  distance  of  the  center  of  gravity  of 
the  aileron  from  the  aileron  hinge  and  represents  the  static  unbalance 
of  the  control  surface.  Positive  values  of  the  parameter  Xp  repesent 

underbalance,  whereas  negative  values  represent  overbalance.  For 

x^  = 0.8  and  for  the  assumed  ratio  of  aileron  mass  to  wing  mass,  values 

of  Xp  ranging  from  -0.01  to  0.01  represent  locations  of  the  center  of 

gravity  of  the  aileron  of  about  20  percent  of  the  aileron  chord  forward 
and  rearward  of  the  hinge.  For  x-j_  = 0.5,  the  same  variation  in  Xp 

would  represent  locations  of  the  center  of  gravity  within  5 percent  of 
the  hinge.  It  should  be  noted  that  in  the  torsion-aileron  case,  a com- 
plete static  balance  in  the  usual  sense  (xp  = 0)  is  not  quite  sufficient 

to  eliminate  mass  coupling.  Actually,  complete  static  balance  against 
rotation  implies  rp2  + 2(x^  - xQ)xp  = 0 so  that  a small  amount  of  over- 
balance (xp  < 0)  is  required.  In  the  present  case,  with  rp2  = 0.004, 
x-l  = 0.8,  and  xQ  = 0.4,  a value  of  Xp  = -0.005  would  be  necessary  to 
balance  the  aileron  against  rotation. 

ftie  parameter  xQ  gives  the  location  of  the  elastic  axis  of  the 

wing,  whereas  ra2  is  the  radius  of  gyration  of  the  wing  referred  to 
the  elastic  axis.  In  the  calculations  for  the  torsion-aileron  case  the 
values  xQ  = 0.4  (elastic  axis  at  40  percent  chord)  and  ra2  = 0.25 
have  been  employed. 


DISCUSSION  OF  RESULTS 


Results  of  the  calculations  are  presented  in  figures  1 to  8.  Fig- 
ures 1 and  2 represent  a general  treatment  with  regard  to  Mach  number. 
Generally,  the  regions  of  flutter  are  the  areas  bounded  by  loops,  either 
open  or  closed.  Where  only  a single  curve  is  shown,  the  flutter  region 
is  the  area  below  the  curve.  For  these  figures  a particular  set  of 
parameters  (x]_  = 0.8,  rp2  = 1/250,  n = 157.08,  % = 0,  and  6p  = o) 
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has  been  selected  and  calculations  performed  for  both  the  bending- 
aileron  and  torsion-aileron  subcases  for  a range  of  Mach  numbers  from  0 
to  2.0.  Figures  1 and  2 present  results  for  the  coupled  flutter  of  the 
wing-aileron  system.  Figure  3 presents  some  additional  results  for  the 
incompressible  case.  Figure  4 presents  results  associated  with  single- 
degree-of -freedom  flutter  of  the  aileron.  For  figures  5 to  8 variations 


have  been  made  in  the  parameters  x^, 


r 2 
0 ' 


p,  g^,  and  gp  and  calcu- 


lations performed  at  isolated  Mach  numbers  to  determine  the  effect  of 
changes  in  these  parameters.  For  convenience  of  reference  to  the  fig- 
ures and  to  show  better  the  range  of  variables  considered  in  the  calcu- 
lations , table  I has  been  prepared. 


For  the  figures  relating  to  the  bending-aileron  case  the  ordinate 
is  a parameter  bo^^cj  the  abscissa,  the  static  unbalance  parameter  Xp 

or  the  frequency  ratio  csp^.  For  the  torsion-aileron  cases  a param- 
eter bo^c  has  been  plotted  against  Xp  or  cup^o^. 

In  figure  1 results  are  given  for  the  bending-aileron  case  and  are 
presented  for  various  values  of  Mach  number,  The  figures  indicate  the 
existence  of  a flutter  region  extending  between  an  upper  and  lower  value 
of  the  parameter  boc^c  (at  M = 0 it  is  expedient  to  employ  bo^v). 

On  the  left  side  of  figure  1,  for  which  the  flutter  parameter  is  plotted 
against  xp  for  the  case  of  zero  restraint  (a>^  = 0) , the  most  striking 

result  is  the  indication  of  a reverse  effect  of  static  mass  balance  at 
the  transonic  speeds.  Whereas  at  M = 2.0  and  also  in  the  incompress- 
ible case,  the  flutter  area  lies  almost  entirely  in  the  region  of  under- 
balance,  as  a Mach  number  of  1.0  is  approached,  the  flutter  region  shifts 
to  indicate  flutter  only  for  the  condition  of  overbalance.  It  may  be 
noticed,  however,  that  the  lower  boundary  of  the  flutter  region  appears 
at  a rather  high  value  of  bo^/c}  this  may  be  a fortunate  circumstance 
in  that  practical  values  of  this  parameter  may  lie  below  the  lower  branch. 
On  the  right  side  of  figure  1,  the  plots  of  the  flutter  parameter  against 
frequency  ratio  ojp  indicate  that  at  a sufficiently  high  value  of 

this  frequency  ratio,  depending  on  the  value  of  the  unbalance  parameter, 
flutter  is  completely  eliminated. 


The  results  given  in  figure  1 relate  only  to  the  coupled  bending- 
aileron  flutter.  At  transonic  speeds  additional  roots  appear  which  are 
associated  with  single-degree-of-freedcm  flutter  of  the  aileron  of  the 
type  discussed  in  reference  13.  These  roots  are  presented  separately 
in  a subsequent  figure. 

In  figure  2 results  are  given  for  the  torsion-aileron  case  and  are 
presented  in  the  form  employed  in  figure  1.  The  results  shown  in  the 
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plots  against  xp  are  generally  comparable  to  those  noted  for  the 
bending-aileron  case  in  that  flutter  is  indicated  at  transonic  Mach  num- 
bers for  the  overbalanced  aileron.  In  contrast  to  the  results  of  the 
bending-aileron  case,  however,  is  the  fact  that  flutter  is  indicated 
for  the  overbalanced  aileron  at  M = 0 as  well.  This  result  for  incom- 
pressible flow  was  further  studied  by  varying  the  mass  ratio  parameter  p. 
Results  are  given  in  figure  3 and  are  discussed  in  a later  paragraph. 

An  attempt  to  study  the  transition  from  M * 0 to  0.7  could  not  be  car- 
ried out  because  of  the  sparse  nature  of  the  available  tables  of  coef- 
ficients. In  the  plots  against  frequency  ratio,  additional  roots  asso- 
ciated with  flutter  of  the  wing  in  pitch  occur  at  transonic  speeds.  It 
might  be  recalled  that  various  studies  have  been  made  for  the  case  of  a 
wing  in  pure  torsion  (see,  for  example,  ref.  14)  and  have  indicated  that 
for  certain  axes  of  rotation  and  certain  Mach  numbers,  flutter  in  this 
degree  of  freedom  can  occur.  Biis  type  of  flutter  was  encountered  in 
the  present  calculations  and  is  indicated  by  the  essentially  linear 
flutter  boundaries  at  Mach  numbers  of  10^9>  10/8,  and  10 jj  (for  the 
latter,  at  high  values  of  only).  At  high  values* of  the  fre- 

quency ratio  these  flutter  boundaries  approach  the  result  of  the  solu- 
tion which  would  be  obtained  by  consideration  of  the  single-degree-of- 
freedom  (in  pitch)  equation. 

As  in  the  bending-aileron  case,  a sufficiently  high  value  of  the 
frequency  ratio  serves  to  eliminate  the  coupled  tor si on -aileron 

flutter  with  the  possible  exception  of  those  cases  which  appear  to  be 
identified  with  single -degree -of -freedom  flutter.  The  required  value 
of  is  rather  high,  however,  and  it  might  be  more  practical  or  more 

desirable  to  eliminate  the  flutter  by  proper  choice  of  bo^/c. 

As  in  figure  1,  the  results  associated  with  single -degree- of - 
freedom  flutter  of  the  aileron  have  not  been  presented  and  are  dis- 
cussed together  with  those  of  the  bending-aileron  case. 

Because  of  the  fact  that,  in  the  tor si on -aileron  flutter  calcula- 
tions at  M = 0,  flutter  was  indicated  for  a highly  overbalanced  aileron, 
this  case  has  been  given  further  consideration.  An  investigation  of  the 
effects  of  various  parameters  indicated  that  the  range  of  values  of  the 
unbalance  parameter  Xp  over  which  flutter  could  occur  for  the  zero- 

frequency  case  was  most  strongly  dependent  on  the  mass  ratio  param- 
eter u.  The  results  of  calculations  for  various  values  of  u are 
shown  in  figure  3*  A decrease  in  u Is  seen  to  result  in  the  division 
of  the  flutter  area  into  two  regions.  As  u is  further  decreased,  one 
of  these  regions  either  vanishes  or  recedes  to  greater  negative  values 
of  Xp  than  were  considered  in  the  calculations.  Flutter  is  still  indi- 
cated, however,  for  a certain  range  of  overbalance. 

The  results  associated  with  single -degree -of -freedom  flutter  of  the 
aileron  which  arose  in  connection  with  both  the  bending-aileron  and 
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torsion-aileron  cases  are  presented  in  figure  4.  The  plots  on  the  left 
show  results  obtained  from  solution  of  the  coupled  bending-aileron  equa- 
tions, while  those  on  the  right  arose  from  the  torsion-aileron  case. 

Note  that  for  the  statically  balanced  aileron  (x^  = 0)  the  two  cases 

yield  identical  results  at  a given  Mach  number,  nils  result  is  the  same 

as  that  obtained  from  solution  of  the  single -degree -of -freedom  equation 

for  flutter  of  the  aileron.  For  values  of  xQ  other  than  zero  the 

I3 

coupled  equations  led  to  slightly  different  curves,  thus  some  coupling 
with  the  wing  motion  is  indicated.  As  mentioned  previously,  the  pure 
single-degree-of-freedom  flutter  of  a control  surface  has  been  treated 
in  reference  13  where  it  is  shown  to  be  highly  dependent  on  altitude, 
moment  of  inertia  of  the  control  surface,  and  structural  damping. 

Figures  5 to  8 give  the  results  of  some  sample  calculations  made 
for  other  values  of  the  control-surface  parameters  and  the  mass  ratio 
parameter  p.  Figure  5 presents  results  for  the  bending -aileron  case 
for  the  values  p = 157*08  and  x-^  = 0.8,  as  treated  in  figure  1, 

with  additional  values  of  radius  of  gyration  parameter  r^2  of  l/l80 
and  lj60.  Consideration  of  figure  5 and  comparison  with  figure  1 indi- 
cates that  an  increase  in  r^2,  corresponding  to  an  increase  in  the 

moment  of  inertia  of  the  aileron,  at  the  higher  Mach  numbers  tends  to 
shift  the  flutter  area  more  into  the  region  of  overbalance,  thus  the 
flutter  area  is  extended  over  a wider  range  of  values  of  the  unbalance 
parameter.  An  increase  in  r g2  also  extends  the  flutter  region  to 

slightly  higher  values  of  ojg j u^.  At  M = 10^9  the  major  effect  of  an 
increase  in  r^2  is  in  narrowing  the  flutter  region  with  respect  to  the 
parameter  ba^jc. 

Figures  6(a)  and  6(b)  show  results  for  the  bending-aileron  and 
torsion-aileron  cases,  respectively,  for  the  set  of  control-surface 
parameters  treated  in  figures  1 and  2 at  a lower  value  of  the  mass  ratio 
parameter  p.  Comparison  of  figures  6(a)  and  6(b)  with  the  cases  for 
lower  mass  ratio  of  figures  1 and  2 indicates  that  at  the  higher  Mach 
numbers  a decrease  in  the  mass  ratio  results  mainly  in  raising  the  values 
of  bo^/c  and  bo^/c  which  bound  the  flutter  region.  At  M = 10/9  a 

decrease  in  the  mass  ratio  shifts  the  flutter  area  farther  into  the 
region  of  overbalance. 

Figure  7 shows  results  at  three  supersonic  Mach  numbers  for  the 
bending -aileron  analysis  of  a large-chord  control  surface  (x^  = 0.5). 

For  this  case  there  is  the  same  reverse  effect  of  mass  balance  at  tran- 
sonic speeds  as  was  noted  for  the  case  where  x^  = 0.8,  and  even  at  the 

higher  Mach  numbers  the  flutter  area  extends  far  into  the  region  of 
overbalance.  With  regard  to  the  plots  against  frequency  ratio,  the 
flutter  region  vanishes  for  values  of  near  unity. 


Figure  8 treats  the  effect  of  structural  damping  for  the  bending- 
aileron  case  and  for  the  set  of  control- surface  parameters  considered 
in  figure  1.  Comparison  of  figure  1 with  figure  8 indicates  beneficial 
effects  of  damping  at  M = 10/7  through  elimination  of  the  flutter 
region  for  the  overbalanced  condition  and  decreasing  the  value  of  ^p/ufo 

required  to  avoid  flutter.  At  M = 1.0  the  addition  of  structural 
damping  shifts  the  flutter  area  slightly  into  the  region  of  underbalance 
but,  for  the  balanced  aileron,  only  a very  low  value  of  a>p/a^  is 

required  to  eliminate  flutter.  At  M = 0.7,  the  amount  of  damping  con- 
sidered eliminates  the  flutter  region  for  the  condition  of  overbalance 
and  greatly  reduces  the  value  of  cup/o^  required  to  avoid  flutter  for 

the  underbalanced  aileron. 


CONCLUDING  REMARKS 


Results  have  been  presented  of  some  sample  calculations  made  for 
the  bending-aileron  and  torsion-aileron  flutter  of  a wing-aileron  system 
of  infinite  aspect  ratio  for  several  values  of  Mach  number  in  supersonic 
flow  and  also  for  some  subsonic  Mach  numbers.  It  has  been  pointed  out 
that  the  results  obtained  are  based  on  the  analysis  of  a highly  idealized 
configuration  and  may  be  considerably  altered  by  such  factors  as  finite 
span,  thickness,  additional  degrees  of  freedom,  and  other  parameters 
which  were  not  considered.  It  is  of  interest,  however,  to  state  the 
effects  of  certain  of  the  parameters  which  seem  to  exert  a major  influ- 
ence since  they  may  serve  to  indicate  trends. 

With  regard  to  mass  balance  of  the  control  surface,  the  results 
indicate  an  apparent  reverse  effect  of  mass  balance  at  the  transonic 
speeds.  In  many  of  the  calculations,  flutter  is  indicated  for  an  over- 
balanced aileron  near  a Mach  number  of  1,  whereas,  at  higher  and  lower 
Mach  numbers,  overbalance  tends  to  eliminate  flutter.  In  practice, 
however,  it  should  be  possible  to  avoid  flutter  for  the  overbalanced 
condition  by  an  appropriate  choice  of  wing  stiffnesses. 

With  regard  to  aileron  frequency,  a high  value  of  the  ratio  of 
aileron  frequency  to  wing  bending  or  torsional  frequency  appears  bene- 
ficial in  eliminating  coupled  wing-aileron  flutter  and,  in  fact,  a 
sufficiently  high  value  will  eliminate  it,  except  in  those  cases  which 
appear  to  be  identified  with  flutter  of  the  wing  in  pure  pitch.  In 
many  cases,  however,  the  value  of  the  frequency  ratio  required  to  avoid 
flutter  appears  to  exceed  unity  and  is  probably  not  easy  to  attain. 

Structural  damping  of  the  wing  and  aileron  appears  to  be  generally 
beneficial  in  reducing  the  amount  of  balance  required  to  avoid  flutter 
or  in  reducing  the  required  value  of  the  ratio  of  aileron  frequency  to 
wing  frequency. 


Near  a Mach  number  of  1 the  possibility  exists  of  flutter  of  the 
aileron  in  a single  degree  of  freedom.  Hie  sequence  of  roots  in  a par- 
ticular calculation  may  be  confused  by  the  appearance  of  roots  associ- 
ated with  this  flutter,  lliis  type  of  flutter  has  been  treated  elsewhere 
and  can  probably  be  eliminated  by  the  proper  choice  of  such  parameters 
as  moment  of  inertia  of  the  aileron. 

Although  certain  parameters  seem  to  give  beneficial  results,  it 
appears  difficult  to  arrive  at  any  single  practical  criterion  for  the 
prevention  of  flutter  involving  control  surfaces  at  supersonic  speeds. 
Further  work  on  the  subject  seems  obviously  needed  since  for  applica- 
tion to  a particular  design  the  effects  of  parameters  not  considered 
will  have  to  be  taken  into  account  before  more  specific  conclusions  can 
be  drawn. 


Langley  Aeronautical  Laboratory, 
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Figure  1.-  Flutter  parameter  bo^y  c or  buc^ fv 
eter  and  frequency  ratio  for  bending-aileron 
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Figure  2.-  Flutter  parameter  tx%yc  or  bo ^/v  against  unbalance  par  am 


eter  and  frequency  ratio  for  tor si on -aileron  case.  |i  * 157»08j 
xQ  * 0.4  j x1  = 0.8;  ra2  - 0.25;  r„2  - 1/250. 
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(a)  Bending -aileron  case. 
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(b)  Tors Ion -aileron  case. 
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ra2  - 0.25. 


Figure  4.-  Flutter  parameter  against  frequency  ratio  for  roots  associated 
with  single -degree -of -freedom  flutter  of  the  aileron. 


Figure  5»-  Flutter  coefficient  tx%/c  against  unbalance  parameter  and 
frequency  ratio  for  bending -aileron  case.  = 157.08$  ■ 0.8$ 

ro2  =»  l/60  and  l/l80. 
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